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THE ANALYSIS OF SOLAR MODELS — NEUTRINOS AND OSCILLATIONS 

Roger K, Ulrich^, Edward J, Rhodes, Jr.^^ 

Steven Toraczyk , Philip J. Dumont', Wendee M. Brunish' 

Introduction 


The theory of stellar structure and evolution is used to 
calculate the properties of a variety of objects from red giants and 
supernova precursors to white dwarfs and neutron stars. Accurate 
tests of the theory in the context of these applications are 
generally not available. The sun as the nearest star provides a 
unique example of a star which can be subjected to a variety of 
precise tests not possible with remote stars. We will concentrate on 
two of these tests — solar neutrinos and solar oscillations — which 
currently indicate that there is something seriously wrong with our 
standard solar model. Although we do not yet know what the source of 
the error is, it is quite possible that the correction of this error 
will require some modification of the results of other applications 
of stellar structure theory. It now seems unlikely that the 
difficulty with the solar neutrino experiment lies in the experiment 
itself. The combination of the difficulty with the solar neutrino 
flux and the difficulty with the solar oscillation frequencies 
suggests that the solar neutrino problem is a failure of stellar 
structure theory rather than a failure of weak interaction theory, 
although this latter possibility cannot yet ba firmly ruled out. 

In addition to the solar neutrinos and solar oscillations, we 
note that two other tests of the sun yield results which are not 
completely understood. First are the abundances of light elements 
which Indicate that the convective envelope extends to 2.8 x 10^ K, 
whereas the standard solar model^ has this temperature equal to 
1.9 X 10^ K. Blake and Schramm^ have considered the possibility of 
convective overshoot and found that this process cannot explain the 
discrepancy without an ad hoc assumption. Second is the shape of the 
solar surface. Apart from phenomena occurring in the visible layers 
of the solar surface, the apparent solar surface corresponds to an 
equipotential surface and provides information concerning the angular 
distribution of matter in the solar interior. The standard solar 
model presumes this distribution to be spherically symmetric. The 
measurements by Dicke and Goldenberg^ have been interpreted by Dicke^ 
to imply a rotating Internal deviation from spherical symmetry. The 
measurements by Hill and Stebbins^ are not in contradiction with 
Dicke* s conclusion, although they are not supportive either. If 
confirmed, Dicke's work would indicate a rather nonstandard property 
of the solar interior. 

The solar neutrino dif f Iculcy^’^ is the best known example of a 
test of the standard solar model which has failed. For a period it 
was common for workers in each of the three areas related to the 
problem — stellar interior theorists, particle physicists and 
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experimental physicists ~ to hope and occasionally believe that the 
solution lay in the other fellow's camp. The experiment has survived 
a very rigorous examination and seems an unlikely place to find the 
cause of the discrepancy. Particle physics, particularly in the area 
of neutrino oscillations, could still play a role in the solution, 
but as pointed out by Bahcall et al.^. oscillations among three 
.^articles would be required and even then some discrepancy would 
remain. Within the context of the standard solar model there are 
uncertainties which do not relate to the assumptions defining the 
Dtandard model which nonetheless lead to uncertainties in the 
calculated neutrino flux. Chief among these uncertainties are the 
possible errors in the calculated opacities, in the measured nuclear 
cross sections and in the abundances of the heavy elements. Bahcall 
e t al . ' » ° have examined the uncertainties and give a predicted 
neutrino counting rate of 7.2 ± 3.3 x 10“^® captures/target atom/s 
(this unit is called the solar neutrino unit or SNU). The 
uncertainty is specifically interpreted as a 3o limit with the 
definition of 3o meaning “if a result falls outside the 3o range, 
then a mistake has been made," when 3o is applied to a theoretical 
result such as the opacity. Since inld-1968, the theoretical 
predictions have fallen within this range, although both the 
prediction and its error have varied. Prior to 1968, the nuclear 
reaction rates were too poorly known to permit an accurate 
prediction. The discrepancy between these predictions and the 
observations (which now yield 2.2 ± 0.3 SNU) has been relatively 
unchanged since the announcement of the earliest experimental results 
in 1968.9 

We wish to emphasize in our discussion the new field of solar 
seismology. This method of studying the solar interior Involves the 
measurement of the structure and frequencies of global solar 
oscillations. The oscillations have been most readily detected 
through the measurement of line-of-sight velocities'^, although there 
have been some measurements of solar limb position variations. The 
largest amplitude oscillations have periods near five minutes. So 
far the oscillations have been studied in three ranges of spatial 
structure — nearly radial, highly nonradlal and most recently 
intermediate size. The elgenmodes are classified according to the 
spherical harmonic T^(6, 9), which describes the angular distribution 
of velocity and according to the distribution and number of nodes in 
the radial component of the velocity eigenfunction. The eigenvalue 
associated with this latter classification corresponds roughly to the 
principal quantum number of atomic physics and ws use the quantity n 
to denote this number. The nearly radial modes correspond to degrees, 
or i values, between 0 and 4, the highly nonradlal modes have i 2 150 
and the intermediate modes have 4 < l < 150. The highly nonradlal 
modes may not be truly global in the sense that I may not be a good 
quantum number as a result of perturbations produced by the solar 
convection zone. The nearly radial oscillations have been detected 
for individual values of n and Jt so that they must be global in 
character. These modes penetrate to the center of the sun and measure 
an average of the sound velocity throughout the solar interior. In 
fact, the spacing of the frequencies from one value of 


n to the next depends on the Inverse of an Integral of dr/c where c 
Is the sound speed. This Integral tends to emphasize the solar 
surface regions, but the deep solar Interior Is also Included. 

Because periods of oscillation can be measured with high accuracy, 
the constraint on the deep solar Interior Imposed by the nearly 
radial oscillation frequencies Is useful In testing both standard and 
nonstandard solar models. 


Observations 


We will concentrate on the nearly radial and Intermediate modes. 
In order to resolve Individual elgenfrequencles , an observing span In 
excess of 24 hours Is required. Such long duration observations can 
be obtained by connecting together sequences from separate days. 

This method was used by Claverle et al.^^ working from the Canary 
Islands to show first that the global oscillations are detectable and 
resolvable. Their observations were Interrupted regularly by the 
diurnal cycle and the lack of continuity of observation complicates 
the analysis of the power spectrum. The diurnal limitations can be 
avoided by going to the geographic south pole during the Austral 
Summer. Informal reports Indicate that a clear spell In the weather 
lasting up to about a week typically occurs at the beginning and end 
of each Austral Summer. During the second of these clear spells In 
1980, Grec, Fossat and Pomerantz^^ obtained the power spectrum shown 
in Figure 1. Tliese remarkable data are clearly worthy of very close 
analysis. 

The sharpness of the power spectriun peaks and their close spacing 
make a detailed examination of the data difficult when they are 
presented In the format of Figure 1. For this reason, it has become 
common practice to use the nearly regular spacing of the frequencies 
to restructure Figure 1 Intc what has been termed an echelle format. 
The scale of the abscissa la stretched and then chopped Into equal 
length Intervals which are then displaced vertically from one 
another. Figure 2 from Grec, Fossat and Pomerantz^^ shows the 
positions of the elgenfrequencles displayed in an echelle format. 

The length of the strips has been chosen to be 136 pHz, the average 
spacing of the elgenfrequencles. Figure 2 compares the frequencies 
from observations by several groups and the overall agreement is 
good. The Identification of the appropriate l value has been made on 
the basis of comparison to theory. The pattern of the eigenfrequency 
spacing Is a more reliable theoretical result than Is the numerical 
value of either the spacing or zero point of the frequencies. 

Comparison of Theory to Observations 

The Identification of the appropriate value of the principal 
eigenvalue n with each power spectrum peak is required before a 
comparison between theory and observation can be carried out. The 
highly nonradlal modes have frequencies which are so widely spaced 
for differing values ^f n that the Identification of n can be made 
easily. By tracing the locus of frequency versus 1 from large i to 
small. It Is then possible to Identify the appropriate value of n for 
the nearly radial modes. Until recently, the Intermediate modes were 
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Figure 1 . The power spectrum of solar oscillations by Grec, Fossat 
and Pomerantz. 


Frequency wHz 

Figure 2 . The positions of central frequencies for the peaks in the 
power spectrum of Figure 1 displayed in an echelle format 
Figure 1 has been divided into 136 uHz intervals which are then dis- 
placed above each other. The different shaped points represent 
results from different groups. 
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not observed and this locus could not be used as a means of 
Identifying n. Now, however, Duvall and Harvey^ ^ have obtained 
frequencies In the Intermediate 1 range and have confirmed the 
Identification made previously on the basis of the closest fitting 
theoretical frequencies. 

Several theoretical Investigations of the theoretical solar 
elgenfrequencles have been made In recent years. The results have 
generally been close to the observed frequencies (within 10 to 
30 yHz) so that the overall picture of the global solar oscillations 
Is reasonably secure. The observed frequencies are defined to within 
1 to 2 liHz for most of the modes, so the possibility arises that the 
standard model Is In disagreement with the observations. Before a 
valid comparison Is possible, however, we need to be sure that theory 
Is reliable within the context of the standard model. Since the 
frequencies are In the range 2000 to 4000 yHz, we need to determine 
the reliability of the calculations at roughly the 0.1 percent level. 
Ulrich and Rhodes^^ have Investigated this question and have 
concluded that the uncertainties In the physics Input Into the 
standard model are not large enough to resolve the disagreement 
between theory and observation. Figure 3 shows comparisons between 
theory and observation for the l ■ 0 to 3 modes, again using the 
echelle format. The discrepancy between theory and observation takes 
the form of both an error In the zero point of the frequency sequence 
and an error in the spacing. The theoretical zero point Is too low 
and the theoretical spacing Is too large. The outer boundary 
condition of the eigenvalue problem as well as some of the details of 
the structure of the solar convection zone just below the photosphere 
are all uncertain and can alter the comparison. However, the spatial 
coincidence of these uncertainties means that they could introduce at 
most one new parameter. Using this parameter, %ie could adjust the 
spacing between the theoretical elgenfrequencles to match the 
observed spacing. If we carry out that adjustment, then we could 
make the zero point smaller and the discrepancy becomes worse. 
Conversely, If we adjust the boundary condition to obtain a larger 
zero point frequency, then the spacing becomes larger and again we 
are no closer to a good fit. The range In the uncertainty In such 
physical parameters as the opacity and nuclear reaction cross section 
Is also not large enough to permit a resolution of the discrepancy. 
Consequently, Ulrich and Rhodes concluded that the discrepancy 
between theoretical and observed frequencies of global solar 
oscillations constitutes a significant failure of the standard solar 
model which should be considered along with the solar neutrino 
problem In the search for a modified solar model. 

Nonstandard Solar Models 


In the course of the search for a solution to the solar neutrino 
problem, several nonstandard solar models have been proposed. 

Table 1 lists a number of the more popular of these models. The 
altered composition models (high Z and low Z) reflect the possibility 
that the solar nebula could have undergone element segregation In the 
condensation phase of evolution and this segregation may not have 





been smoothed out by a fully convective phase* Neither supposition 
is supported by our present theories of star formation* but the early 
stages of stellar evolution are sufficiently complex that such 
segregation cannot be ruled out* Although as we will see below both 
the high and low Z models fail* we feel the idea of element 
segregation is attractive* Another version of modified element 
distribution was suggested by Schatzman and Maeder^^ and involves 
diffusive mixing* Their hypothesis is for a small scale turbulence 
field to enhance the diffusion coefficient enough that the gradient 
of hydrogen abundance is smoothed out* This model is attractive 
because of its low neutrino flux; however* there is no known 
mechanism which could enhance the diffusion coefficient enough to be 
consistent with the model* Another version of element segregation 
Involves modifying the Initial distribution of the hydrogen and 
helluTj throughout the sun* This could have occurred if, for example* 
very thick molecular hydrogen mantles were formed on the 
surfaces of grains in the coldest part of the protosolar cloud* 
Segregation between the grains and the residual gas is possible if 
the grains are charged and bound to a magnetic field while the gas is 
neutral* We might expect the solar core to be relatively depleted in 
both hydrogen and heavy elements according to this scenario* 

Finally* wc have considered a truncated Cowllng^^ magnetic field of 
high strength* The central field is large enough that the magnetic 
pressure |Bp/8ir is a few percent of the gas pressure* This requires 
B - 3 X 10® gauss at the solar center* Various authors^' have argued 
that 1) a field larger than 10^ gauss will produce an unacceptably 
large solar oblateness and 2) magnetic buoyancy will cause a field 
larger than 10® to 10^ gauss to rise to the solar surface* We do not 
know of an argument which could counter the first point* although it 
is possible in principle that some combination if dipole and toroidal 
fields could avoid producing significant distortion on the solar 
surface* We believe that a large enough scale field may not be 
subject to the buoyancy considerations* at least for a long enough 
time that nuclear burning can enhance the helium abundance to the 
point where it can counterbalance the magnetic buoyancy* Several 
additional nonstandard models are also listed in Table 1* which 
differ only slightly from the standard model* These are Included as 
tests of the effect of uncertainties in the basic physical Input on 
various derived model properties* 

Tables 2 and 3 give some of the Important derived properties of 
the models listed in Table 1* The neutrino fluxes in the third and 
fourth columns are larger for the standard model than have given by 
Bahcall ^ al*^ because the code used to calculate solar models for 
the purpose of studying elgenf requencles uses a less accurate 
treatment of the abundances of the minor nuclear constituents than 
does the code used for the solar neutrino problem* Ratios between 
the fluxes given in Table 2 for different models should be reliable 
even though the absolute value is high* The frequencies in Table 3 
differ slightly from those plotted in Figure 3 for the standard 
model* This is because Table 3 Includes the results of a new 
radiative Interaction theory* Also* the radius is corrected to the 
more precise value given below* 
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Table 1 


Model Characteristics 


Model Name 

Model Description 

Standard 

Has standard physics, Initially chemically 
homogeneous, no mixing, no magnetic fields. 

High Z 

Has Its Internal heavy element abundance Z equal 
to 1.19 times Zo^ (2^ refers to abundances from 
Ross and Aller^" and is 0.018). The surface 
abundances are normal (Z > Z^) . 

Low Z 

Has Internal heavy element abundances Z ■ 0.3 
Zra« The outer abundances are normal except for 
a small excess In hydrogen relative to helium to 
maintain a mean molecular weight gradient which 
Is normal. 

Low X, Z - 0.01 

The hydrogen mass fraction X Is dropped by 0.04 
and the heavy element abundance Is dropped to 
0.01 Interior to M ■ 0.7 1^. The helium 
abundance Is Increased In the Interior to 
compensate. 

Diffusive Mixed 

Enhanced diffusion mixes fresh hydrogen Into the 
solar core as suggested by Schatzman and 
Maeder. 

Magnetic A 

Has a Cowling magnetic field with Bg ■ 3.16 x 10® 
gauss. The first zero In B Is at r » 0.7 R 0 . 

The field Is assumed to remain zero exterior to 
this point. 

Magnetic B 

Like Magnetic A but with the zero at r ■ 0.5 

Inner 0.05 

The Inner 0.05 of the sun's mass Is artificially 

mixed 

mixed on a time scale rapid enough to homogenize 
the -^He as well as the hydrogen. 

S 34 “ 0.3 KeV-Barnes 

The lower value of the ®He + ^He ^Be + y 
reaction suggested by the Munster group Is 
adopted. 

Enhanced npacity 

The opacity Is artificially multiplied by 1.2 
throughout the model. 

No scattering 

The second vlrlal coefficient due to scattering 

states 

states Is set to zero. 


Model Derived Properties 
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The behavior of Che physical parameters at the base of the 
convective envelope for the different models la not Intuitively 
obvious. This is because Che envelope Is adjusted to force the solar 
model radius to be 6.9625 x cm^^ at an optical depth of for 

normal Incidence rays. Some of the quantities varied, such as the 
opacity, would cause the envelope depth to change If all other 
quantities were held constant. For example, the higher X In the 
envelope of the low X and low Z model would tend to Increase Che 
temperature and decrease the density at a fixed depth because of the 
lower mean molecular weight. However, the total volume available to 
the model Is fixed so that the convective efficiency Is Increased In 
order to achieve a satisfactory model. It Is the Interplay between 
convective efficiency, average hydrogen abundance and the Introduced 
model variations which causes the results to be difficult to deduce 
without the detailed calculations. 

We see from Table 3 Chat the model with a discontinuity In 
hydrogen abundance as well as In heavy element abundance yields 
frequencies in good agreement with observation. Although the general 
pattern of abundance modification which this model represents Is 
attractive on the basis of this result, we believe that It almost 
certainly la not unique. There are two parameters which are 
adjustable — the Initial X abundance gradient and the position 
within the sun where It occurs. Our Incomplete results for other 
models not listed In Tables 2 and 3 Indicate that these two 
parameters can be adjusted simultaneously so that the frequencies of 
the nearly radial modes are unchanged. We note, however, that this 
form of model Is the only one we have yet found which appears capable 
of yielding frequencies In agreement with the observations. 

Finally, wa wish to emphasize the potential Importance of the 
Intermediate degree modes. For a fixed frequency, these have very 
similar eigenfunction structure In the outer parts of the sun, but as 
£ decreases their Inner turning point penetrates progressively more 
deeply. Thus, the comparison of solar frequencies with the 
theoretical models at £ valuea between 0 and 100 can be used to probe 
solar structure In a selective fashion. The frequencies are 
dependent on the sound velocity as a function of depth. Figures 4 
and 5 llluatrate the similarity of the frequency changes as a 
function of £ to the sound velocity changes as a function of r/R^. 
Clearly, there Is great potential for using solar seismology to 
choose among the various nonstandard solar models. 
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